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Model independent reconstructions of dark energy have received some attention. The approach 
that addresses the reconstruction of the dimensionless coordinate distance and its two first deriva- 
tives using a polynomial fit in different redshift windows is well developed [4H6|. In this work we offer 
new insights into the problem by focusing on two types of observational probes: SNela and GRBs. 
Our results allow to highlight some of the intrinsic weaknesses of the method. One of the directions 
we follow is to consider updated observational samples. Our results indicate than conclusions on 
the main dark energy features as drawn from this method are intimately related to the features 
of the samples themselves (which are not quite ideal). This is particularly true of GRBs, which 
manifest themselves as poor performers in this context. In contrast to original works, we conclude 
they cannot be used for cosmological purposes, and the state of the art does not allow to regard 
them on the same quality basis as SNela. The next direction we contribute to is the question of how 
the adjusting of some parameters (window width, overlap, selection criteria) affect the results. We 
find again there is a considerable sensitivity to these features. Then, we try to establish what is the 
current redshift range for which one can make solid predictions on dark energy evolution. Finally, 
we strengthen the former view that this model is modest in the sense it provides only a picture of 
the global trend. But, on the other hand, we believe it offers an interesting complement to other 
approaches given that it works on minimal assumptions. 
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I. INTRODUCTION 

It is very timely to continue to explore the accelera- 
tion of the Universe [I] , as its discovery has very recently 
finally gathered the top recognition of the community. 
The question of whether such evolution is produced by 
a dark energy component or rather a manifestation of 
bold unknown Physics remains one of the most puzzling 
topics in Cosmology Assuming the dark energy sce- 
nario, and beyond theoretical questioning about its phys- 
ical nature, it is most vital to learn to extract the best 
and most relevant information from current astronom- 
ical data. Many methods are strongly tied to the as- 
sumption of a background cosmological model or, equiv- 
alently, of an analytical dark energy equation of state 
(EoS). Such is precisely the case of the reconstruction 
of the Hubble diagram of type la Supernovae (SNela), 
which provided the first evidence for the acceleration of 
the Universe [1] . The diagram presents the corresponding 
distance modulus at different redhsifts and is computed 
from the EoS itself after a convenient integration in red- 
shift. Variations of this procedure apply as well when 
considering other observational tests such as gamma ray 
bursts (GRBs), baryon acoustic oscillations (BAO) or the 
cosmic microwave background (CMB). But even though 
biases on final results may depend on the previously de- 
scribed assumptions, that route, and in particular giving 
a specific prescription of the EoS is still the most popular 
data analysis method, as opposed to model independent 
methods (EoS free prescriptions) for the reconstruction 
of dark energy. Fortunately, literature sources about this 
alternative way to explore the accelerated universe have 



mounted considerably, and thus the state of the art has 
turned promising [3]. 

In [IHS] the authors present a different approach to this 
model independent reconstruction: the dimensionless co- 
ordinate distance and its first and second order deriva- 
tives with respect to redshift are reconstructed using a 
polynomial fit inside an appositely designed redshift win- 
dow. The only underlying assumption is a description of 
the Universe by a generic Friedmann-Robertson- Walker 
(FRW) metric. Moreover, the method does not depend 
on the physical properties of the components of the Uni- 
verse or on the underlying gravity theory. The latter is 
an aspect to be considered only at a second stage, when 
one attempts at the derivation of the behavior of more 
interesting and observable cosmological quantities: the 
dimensionless expansion rate E{z) = H{z)/Hq, the uni- 
verse acceleration (through the well known deceleration 
parameter q(z) = —da/cP), or finally the dark energy 
EoS parameter w{z) = Pde/ Pde- As usual, here a is the 
scale factor, the dots denote derivatives with respect to 
cosmological time, and pde and pde are respectively the 
dark energy pressure and density. 

In this work we aim at a thorough and novel study of 
a number of main aspects of this problem. Firstly, ob- 
servational samples are updated, and results with are 
compared, as far as possible. This lets us verify that, op- 
posite to what one could expect in principle, conclusions 
in general are strongly related to the way observational 
samples have been built on. This is not, thus, a caveat 
of the reconstruction method, and could rather be useful 
at a first observational level so as to discern the require- 
ments a sample needs to be useful for the inference of 
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information about dark energy. Secondly, the validity of 
GRBs in the context of this method is checked. It has 
been pointed out [3] that these data can be a promising 
tool to study cosmology at high redshifts. In this regard, 
our results seem to hint that GRBs are not as compelling 
as desirable, and probably the reasons for this arc intrin- 
sic to the nature of actual GRBs observables. This is in 
agreement with the critical view in [7^ . 

After this, it is argued that the analysis in the model 
independent spirit of [3] requires adjusting some param- 
eters. But typically, this causes noise and large errors 
in some regimes (this seems to be intrinsically related to 
the procedure itself). Thus, the following questions arise 
quite naturally: how can the real cosmological informa- 
tion be read off from these noisy patterns? And how 
and when can the goodness and closeness to reality of 
the reconstructed profiles be ascertained? 

At the next stage, we move on to explore the range of 
application of the method and compare it with the results 
obtained using a specific dark energy EoS. At the same 
time, we investigate how physically relevant parameters 
depend on the errors and noise of data. 

Finally, we conclude that this method just provides an 
idea of the global trend of the parameters. The large 
correlation among them does not allow to obtain inde- 
pendent measures. But, on the other hand, it is admissi- 
ble to report values associated with independent redshift 
bins. 

In Sec. In] we revise the method in its main aspects; in 



Sec. Ill we briefly describe the observational tools chosen 
for the analysis, namely the SNela and GRBs samples 
with their properties; in Sec. |IV| we discuss the main 
results both for real and mock data; in Sec.[V]we present 
results from the binning approach; and in Sec. VI we 
summarize the most important results. 



II. METHOD 
A. Basic quantities 

The use of SNela and GRBs makes us focus mainly on 
the distance modulus /i, which is an observable quantity 
defined as 



^(z) = 51ogio[i?L(z)] +^0 



(1) 



Here is the luminosity distance (see below), the nui- 
sance parameter /io encodes the absolute magnitude M 
of any SNela and the Hubble constant Hq (the current 
value of the expansion rate). For convenience, the speed 
of light has been set to one, and this choice will apply to 
the whole paper. 

In full rigour, though, a proper model-independent 
analysis of the GRBs data cannot be based on the GRBs 
distance modulus, which is not directly observable and 
depends on the cosmology assumed [H]. As we will dis- 
cuss in IIIB we have used a GRBs sample whose cali- 



bration is independent of the cosmological model and is 
based on the use of SNela as calibration tools. 
The luminosity distance can be written as 



Dr 



H^\l + z)y{z). 



(2) 



where the dimensionless coordinate distance y{z) reads 
explicitly 



y{z) = 



dz' 
W) 



(3) 



The determination of the coordinate distances from the 
underlying observational data typically proceeds under 
the assumption of a cosmological model, which reflects 
in the analytical expression of E{z), the dimensionless 
expansion rate. 

The method developed in [4] does not require any as- 
sumption on the cosmological model, or on the expres- 
sion of the dark energy functional. It is possible though 
to determine y{z) in a model independent way. Then, 
as determining the Hubble function H(z), and deceler- 
ation parameter q{z) only involves the computation of 
the derivatives of y{z), those two quantities of interest 
inherit the model independence. 

The method was developed by [4^ and follows a series 
of steps which we describe in the next subsection. 



B. Inputs and fits 

Step 1. The input data needed are the collection of 
{zi , yi , ct|), quantities, where the index i ranges from 
1 io N (the number of objects in the chosen dataset). 
Let us describe those quantities one by one. Clearly, the 
Zi are the redshifts for the dataset objects. Then, the 
yi = y{zi) denote the coordinate distance data values ob- 
tained from the observational quantities (as the observed 
variable is ^i{z), Eqs. ([T]) - ([2| have to be combined to 
give the yi). Finally, the af represent the coordinate dis- 
tance errors. Contributions to this last quantities will 
come both from errors on /i(zj) and from errors on the 
redshifts z;. The sample we have used was not presented 
with errors on the redshifts, but using old dataset val- 
ues, we have verified that the redshift contribution from 
those errors is always negligible when compared to the 
contribution from the errors in the fi{zi). 

Step 2. We define a redshift grid with the values Zj, 
and at each of them we evaluate the coordinate distances, 
their derivatives and the related cosmological quantities. 
The properties of this grid are completely arbitrary, but 
we have checked this does not affect final results; in our 
case the grid will range from until the maximum red- 
shift z^ax reached by the chosen dataset, with an homo- 
geneous spacing of 0.005. 

Step 3. Finally, we have to perform a weighted poly- 
nomial fit of {zi , yi), with the inverse errors l/crf as sta- 
tistical weights. The fit is not, though, a global one. We 
rather do a number of individual fits in separate redshift 
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windows defined by a length Af and a center zq. For 
each fit just those data points falling in each window are 
taken, and the central redshift moves iteratively along 
the grid values Zj. At both ends of each window, we 
apply in addition a Gaussian tapered region with disper- 
sion cTg and extending out to 2a g. This makes all the 
points in this region have their weights reduced by the 
value of the Gaussian wing at that point. This way, fluc- 
tuations and discontinuities are controlled when moving 
from one window to another. Then the fit is performed 
after changing the variables from zi to = Zj ~ zq, so 
that the coordinate distance values and its derivatives at 
zq = Zj are simply the values of the coefficients of the 
polynomial evaluated in a; = 0. For example, when using 
a second order polynomial: 



y{x) = flQ + aix + a2X^, 



we will have 



2/(0) = ao, 2;'(0)=ai, y"(0) 



02; 



(4) 



(5) 



where the prime denotes differentiation with respect to 
X (or, equivalently, to redshift, in this case). All the cos- 
mological quantities depend on combinations of the re- 
spective (y(0), y'(0), y"(0)), and the errors can be easily 
evaluated by straightforward application of error propa- 
gation theory. 



C. Outputs and derived quantities 

The method described in the previous subsection pro- 
duces a set of j polynomial-reconstructed dimensionless 
coordinate distance values with their first and second or- 
der derivatives values (here j is an integer depending on 
the redshift grid range and spacing). This allows giving 
the next step in the reconstruction. 

Assuming an FRW metric and no contribution from 
spatial curvature, one can easily obtain estimations for 
more cosmologically interesting quantities. One of them 
is the expansion rate, which from Eq. ([3| is easily deriv- 
able to be 



E{z) = (y')- 



(6) 



where the prime indicates differentiation with respect to 
redshift. Likewise, the deceleration is given by 



q{z) = -^ = -l-{l + z)^- 



(7) 



This expression is obviously valid for any cosmological 
scenario based on an homogeneous and isotropic expan- 
sion. If we now also assume a theory of gravity (in this 
case General Relativity), we can use the Friedmann equa- 
tions: 



a 
a 



47rG 



(Pm + PDE + iPDs) , 



(8) 
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(Pm + PDe) : 



(9) 



where Pm is the mass-energy density of non-relativistic 
matter, pde '^s. the mass-energy density of the dark en- 
ergy, and pde is the pressure of dark energy. It is 
straightforward to derive from them the expression of 
the dark energy equation of state w. Using the standard 
definition of the critical density at the present epoch, 
po^c — 3i?Q/(87rG), the dark energy pressure can be 
shown to be 



Pde(z) = po,, 



Po,c 



E\z) 



-iy' 



{2q{z) - 1) 



2 . ^ y" 



; (10) 



while the dark energy density is 

Pde{z)^ Po,c[iy'r^ ~n„,il + zf] , (11) 

where Qm — Pm/po,c is the fractional contribution of 
non-relativistic matter to the critical density at zero red- 
shift, and the term (1 -|- z)^ defines its time evolution. 
Then, as the dark energy EoS parameter w{z) is de- 
fined as the ratio of the dark energy pressure to density, 
Pde{z) / P£)e{z), it finally results in the expression 



w{z) 



-[l + (2/3)(l + z)y"(yr 
l-a„(l-t-z)3(y')2 



(12) 



Note that is clearly involved in the determination of 
the w{z) parameter and of its related error; here we use 
the value corresponding to the wcdm+sz+lens case from 
WMAP7-year [9j, which has r2,„ = 0.259ij]:[!^^. Nev- 
ertheless, this particular choice does not have a strong 
impact in the conclusions to be reached. 

To close this section let us mention that errors on all 
the previous quantities are of course derived from errors 
from y, y' and y" applying the usual error propagation 
rules. 



III. CURRENT DATA 

A. Supernovae 

Our main reconstruction tool is one of the most recent 
type la supernovae sample available, the Union2 sam- 
ple described in [10]. This compilation is the result of 
a new dataset of low-redshift nearby-Hubble-flow SNela 
and is built with new analysis procedures suitable for 
working with several heterogeneous SNela compilations. 
It includes the Union data set from [llj with six added 
SNela first presented in [10^ , along with SNela from [H] , 
the low-z and the intermediate-z data from [T3] and [13] 
respectively. After the application of various selection 
cuts to create a homogeneous and high signal-to-noise 
data set, A/sNcia = 557 SNela events distributed over the 
redshift interval 0.015 < z < 1.4 were obtained. 
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The use of the Union2 SNela data rests on the defini- 
tion of the distance modulus as given above. Given the 
heterogeneous origin of this dataset, and the procedures 
described in 11] for reducing data, it turns out we need 
to fix the value of /ip for converting the modulus distance 
/i(z,/io) to the coordinate distance y{z). One way would 
be using just low redshift SNela to find the best fit value 
for Hq starting from the approximate expression valid for 
this range, Dl ~ z/Hq; alternatively, one could assume 
an independent measurement of the same quantity, as in 

m- 

Nevertheless, we have chosen a different way to pro- 
ceed: as we will try to compare results from SNela with 
results from GRBs, choosing the same unique value for 
Hq makes sense, even though it is not fundamental for 
the final results. For this reason, as the GRBs calibration 
has been done assuming a well fixed value of Hq = 70 
(km/s)/Mpc [HI [TS], we have decided to use the same 
value of Hq for the SNela (which at the same time we 
have tested to give a very good fit to the Hubble dia- 
gram of SNela at low redshifts). 



B. Gamma Ray Bursts 

Working out and interpreting results from a GRBs 
analysis is not an easy task, as the errors on their ob- 
servable quantities are much larger than those for SNela, 
besides the fact their source mechanism is still not well 
understood. For this reason, choosing a good GRBs sam- 
ple is crucial; we have chosen to work with the sample 
described in |17) . There the authors performed a new cal- 
ibration procedure on the widely used GRBs sample from 
[15], which matches perfectly some of the requirements 
in this paper. The possibility of a comparison between 
SNela and GRBs is strictly related to the building of 
a Hubble diagram for GRBs too, which is an extremely 
difficult task because GRBs are not standard candles, un- 
like SNela. To create an Hubble diagram for GRBs, one 
has to look for a correlation between a distance depen- 
dent quantity and a directly observable property. Start- 
ing from some of the many correlations that have been 
suggested not long ago, in [TS] the authors created a sam- 
ple of 69 GRBs in the redshift range 0.17 < z < 6.6 for 
which the Hubble diagram is well settled. 

Then, in |17) that sample was updated in many as- 
pects, the main one being the test of a new method for 
the calibration of GRBs based on the non assumption of 
a priori cosmological model. Such a model independent 
calibration is built on the idea that SNela and GRBs at 
the same redshift should exhibit the same distance modu- 
lus. In this way, interpolating the SNela Hubble diagram 
gives the value of ij.{z) for a sub-sample of GRBs which 
lies in the same redshift range. This sub-sample can be fi- 
nally used to calibrate the well known GRBs correlations 
and, assuming that this calibration is redshift indepen- 
dent, it can be extended to high redshift GRBs. With 
this procedure, in TT the authors were able to convert 



the [T5] sample to a new one, with the same number of 
objects but with SNela-calibrated GRBs distance modu- 
luss. 



IV. RESULTS 

We remind that the main objective of this kind of ap- 
proach is not to obtain independent measurements of co- 
ordinate distances and their related quantities, but rather 
to find out the global trends underlying the data. More 
exactly, we aim to specify a possible way of proceeding 
for optimizing the use of this approach. In particular, 
we will examine the effects produced by changes in the 
fitting windows which should enable to detect and dis- 
tinguish global and local trends. 

On the one hand, as it is well known, using larger fit- 
ting windows leads to more robust fits, as more points are 
included in the analysis; but, at the same time, a strong 
correlation arises among the various values one obtains, 
as the fitting windows clearly overlap. 

On the other hand, using narrower fitting windows, 
makes it possible to find local (in a redshift sense) prop- 
erties of the global trend, and with more redshift resolu- 
tion. The negative counterpart is that smaller windows 
also imply working with fewer data points, and this ac- 
centuates noisy features in the estimated parameters. 

If one assumes implicitly the physically reasonable hy- 
pothesis that the coordinate distance y{z) is a smooth 
function of redshift (as one does in this kind of approach) , 
then one can compare large (global) and small (local) fit- 
ting windows results. In that way, one can explore the 
global trend and the characteristics that depend on local 
irregularities coming from, for example, mixing hetero- 
geneous sub-samples, even if they are treated appropri- 
ately, as is the case of the Union2 sample, or containing 
intrinsic and physically important information about the 
universe dynamics. By comparing results from large and 
small fitting windows, one can detect more clearly recon- 
struction features due to the limitations of the method 
rather than to an intrinsical relation to the physics dy- 
namics. In this work we study which is the best range 
to obtain this information, and we also gather some fur- 
ther understanding about the limits of the method pre- 
sented: we basically conclude the limits depend on the 
characteristics of the datasets and the specifications of 
the fitting windows. We have also performed an indepen- 
dent binning sketch, which has consisted in dividing the 
data points into a number of bins which do not overlap, 
and obtaining independent estimations of the required 
parameters afterwards. 



A. SNela: comparison with literature 

To avoid misunderstandings, let it be clear that there 
are a number of facts that preclude a very close com- 
parison of our results with those in To begin with. 



5 



our SNela sample has two main differences with respect 
to the samples used in [6]. First of all, the sample we 
consider is made up of 557 objects as opposed to the 192 
in Davis' dataset. This mainly implies a higher mean 
number density for the data points which should trans- 
late into a better knowledge of the global tendency of the 
parameters under study. This newer sample displays a 
quite larger number density for z < 1, but no consider- 
able changes occur for higher redshifts values, so one still 
has to be careful about conclusions based on the results in 
this range. Second, the sample we use is made of different 
heterogeneous SNela sub-samples and, even though data 
reduction has reduced differences considerably, it is pos- 
sible that those disparities influence the local properties 
of our derived quantities. 

A quite important new point in our analysis is that we 
have noticed that some remarkable oscillations that arise 
at some levels of the reconstruction are not directly due to 
the fitting algorithm, but mainly related to the presence 
of SNela with a high relative error. To back up this 
view, we produced a cut version of the total Union2 data 
sample by deleting all the data points whose error was 
larger than < err > -f 2 Gem where < err >— 0.225 is 
the mean of the data error distribution and (Terr = 0.133 
is its standard deviation. With this procedure we cut 30 
points (all having in common a high relative error) and 
the oscillatory (noisy) features disappeared. 

We performed different runs changing both the length 
of the fitting windows. A/, and the dispersion ag of the 
Gaussian tapered region needed for calculating the fit 
weights. We regarded as good results just those obtained 
when the number of data-points inside any fitting win- 
dow and out of the Gaussian tapered region (i.e. with 
unchanged error weight) was > 10, so to assure a minimal 
sufficient number of data points as compared to the num- 
ber of fitting parameters. Our choice is slightly stricter 
than the one proposed in [3]. 

Fig. [ija) shows a reconstruction of y derived from a 
fit with Af = 0.6 and a dispersion ag = 0.02 of the 
Gaussian tapered region (i.e. the same values used in [5] 
and regarded there as the best ones) . The trend is clearly 
very smooth and quite similar to the one in [B] ; all in all, 
it can be considered as an admirably good reconstruction. 
We detect, though, that for redshifts z > 0.8 there is a 
slight difference between the reconstructed dimensionless 
coordinate distance and the theoretically expected one 
(derived using some updated cosmological models from 
[5] with respect to the one used in [S]), with a systematic 
lower trend at high redshifts. Given the location (redshift 
range) of the deviation, we believe it can put be down to 
the low number density of SNela points in this region. 

The reconstruction of the first order derivative of di- 
mensionless coordinate distance y' is clearly less satis- 
factory. We plot it in Fig. [l](b) using the same redshift 
range as for the corresponding figure in [6] . Both studies 
display in a broad sense similar irregularities for z > 0.6, 
but a rise at high redshifts occurs in our case; it is either 
a hint of a possible systematics in the way high redshift 
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Figure 1. Comparison of y{z), y'{z) and y"{z) with A/ = 0.6 
and ag — 0.02 from a polynomial analysis with a WMAP7 
ACDM (STi™ = 0.266, solid line) and with a WMAP7 {+ 
WMAP7.2) quiessence model {Qm = 0.27, w = -1.1, dashed 
line). 



data have been obtained, or it is rather a natural conse- 
quence of the low number of data points per bin used in 
the interpolation, or just a mixture of those such prob- 
lems. No way of disentangle the two effects is possible at 
this level. 

If we compare this trend in y' with the results about 
y, it follows immediately that the redshift range where 
y' rises is the same where y deviates more from the the- 
oretical expectations. We underline that the theoreti- 
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cal expectations (plotted in the same figures) are only 
considered as a possible comparison term and a tool for 
detecting possible local or global properties of the under- 
lying trend. Their use does not mean they are considered 
important as a fiducial cosmological model. 

For what concerns the second order derivative, y", 
Fig. [ijc), we clearly observe a very noisy behavior at 
redshift z > 0.4. This quantity seems more sensitive to 
oscillations in the number density in SNela sample. We 
reach this conclusion because while [6] this trend looks 
quite homogeneous, in our reconstruction it is more ac- 
centuated so we are tempted to put it down to the hetero- 
geneous nature of the Union2 sample. Thus, even though 
we are working with a larger SNela sample, error bars 
are quite comparable and, in our case, even larger at 
high redshift. So, the first lesson one can derive is that 
denser and larger (in a number data point sense) samples 
might not improve the use of the model independent re- 
construction method we are analyzing. This can be seen 



as a limitation of the method itself, and in IV D we look 
deeper into this problem. 



B. SNela: general overview 

Continuing with our analysis, we have also considered 
differences arising from the use of different window width 
values. After several preliminary tests, we have eventu- 
ally considered as representative results from A/ = 0.6 
(which is the width used in [B]), and from another two 
widths: A/ = 0.4 and A/ = 1.0. 



1. How small a window? 

As we are working with a data set that contains a 
statistically significant number of objects, one may be 
tempted to try and improve the analysis by using smaller 
redshift windows for the fit, as, in principle, one would 
still have a sufficient number of points. This is, however, 
not true mainly due to the non uniformity of redshift dis- 
tribution of our SNela sample. We have a large number 
of objects at small redshifts but too few points are at 
high redshifts, and this makes the noisy trends and the 
errors grow in this range. It is thus not surprising that 
more stringent constraints on our derived quantities at 
large redshifts cannot be obtained. 

Moreover, when the fitting window width becomes 
smaller, the extremal intervals (i.e. intervals at very low 
and high redshift) become precisely the most problematic 
cases. We can see this in Fig.l2je) by paying attention to 
the case of the reconstructed y^^. Note that the errors get 
minimal in the mid redshift range, but then grow again 
for low (large number of data points) and high (small 
number of data points) redshifts. On broad terms, at 
low redshifts the number of datapoints in typical win- 
dows is high, and so we have small noisy oscillations and 
biased estimations of the reconstructed parameters, but 



estimations also depend on the Gaussian tapered region 
dispersion, as- In contrast at high redshifts, few points 
will be found in our window, and constraints are poor: 
many oscillations and large error bars occur. 



2. How large a window? 

On the other hand, taking into account that our main 
goal is to guess the global trend of our parameters, 
we considered a somehow completely opposite approach. 
Specifically, we studied the effect of enlarging the fitting 
window as much as possible so as to detect a global and 
unified trend of the data (even though this would cause 
overlaps in the fitting windows and a related growth in 
the correlation). We detected that large fitting windows 
give very smooth trends, mainly because the data are, 
in some measure averaged. But at the same time these 
results are strongly biased (as already reported by 0), 
because the polynomial fit describes very poorly the gen- 
eral trend of our quantities. More specifically, a second or 
third order polynomial can fit very well the dimension- 
less coordinate distance in the redshift range depicted 
by data, but it would give strongly unphysical trends 
when extended to higher (out of sample) redshifts as, for 
example, a negative coordinate distance. For this rea- 
son, we do not present here cases with widths wider that 
Af — 1.0, because this is the largest amplitude for which 
such a bias has a negligible effect on our reconstructed 
quantities. 



3. Parameter estimates 

Let US review now our estimations of several parame- 
ters, obtained under the considerations described above. 

a. Hubble parameter The reconstructed Hubble di- 
agram hardly changes for the three different windows 
lengths considered. As we show in the zoomed vision 
of Fig. [2](b) , differences are negligible and no significant 
discrimination is possible when compared with Union2. 
We have to recall the good behaviour of the reconstruc- 
tion of this quantity. The reconstructed Hubble param- 
eter. Fig. [2](d), is practically independent of the fitting 
window for z < 0.5, and in this redshift range it fits 
quite well the observational data extracted from differ- 
ential ages of passively evolving galaxies and discussed in 
[ig. The Hubble constant (i.e. Hq = H{z = 0)) is found 
to be: 



70.63 ± 1.85 for A/ = 0.4, 
Ha = { 71.14 ± 1.34 for A/ = 0.6, 
71.04 ± 1.02 for Af = 1.0 . 



(13) 



in units of (km/s)/Mpc. Clearly, differences among the 
three fitting windows are very small, just notice that the 
larger the window, the lower the error on Hq. We can also 
see that the three trends are quite similar up to 2 « 0.5, 
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Figure 2. Second order polynomial reconstruction with SNela (Union2). Coordinate distance 3/(2) and its derivatives {y'{z) and 
y"{z)), and the derived cosmological quantities; distance modulus ^j.{z), Hubble function H{z), deceleration parameter q(z) 
and the dark energy EoS w{z). Different colors are used for different fitting window A/: from lightest to darkest for A/ = 0.4, 
A/ = 0.6 and A/ = 1.0 respectively. 
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while significant deviations appear for to high redshift 
data. 

b. Deceleration parameter The deceleration param- 
eter q{z), Fig. [2jf) is influenced by the noisy features in 
y"{z), which become rather important at high redshifts. 
Differences among the three fitting windows considered 
become manifest when one calculates the present value 
of q{z) and zt (transition redshift from deceleration to 
acceleratiorQ Our estimates are: 

{-0.477 ±0.391, zt-0.31 for A/ =0.4, 
-0.652 ± 0.185, Zt - 0.38 for Af = 0.6, (14) 
-0.645 ± 0.091, zt-0.51 for A/ = 1.0. 

Just for comparison, the WMAP ACDM model plotted 
in Fig. [3] has qo — —0.601 and zt — 0.77, while the 
WMAP quiessence model q^ — —0.704 and = 0.74. 
In this sense, the reconstruction at low redshifts seems 
to be consistent with such cosmological models, so we 
have another proof of the consistency of the method in 
that regime. 

c. EoS parameter It is not an easy task to argue 
whether the EoS parameter w{z) has a constant or a dy- 
namical behaviour (as it is possible to verify by having a 
look at Fig. |4| : for high redshifts the noise and the errors 
make it impossible to draw any strong conclusion. Nev- 
ertheless, we can offer estimates for the present values: 

f -0.874 ± 0.369 for A/ =0.4, 
Wo -1-025 ±0.211 for A/ = 0.6, 
[-1.019 ±0.154 for A/ = 1.0. 

Note that those are all compatible with both the WMAP 
ACDM and quiessence model. Another interesting value 
is ^0.5 = w{z = 0.5): following [TS] it is worth con- 
sidering an alternative parametrization for the dark en- 
ergy EoS, as opposed to the most used [20] Chevalier- 
Polarski-Linder (CPL) modeQ [211 EZ]. This scenario 
(specifically, a reparametrization of the former) addresses 
the values wq and W0.5, instead of the asymptotic value 
Wa = w{z — >■ 00) (see footnote). This pair of values has 
the advantage of lower correlation and informs us about a 
redshift range well accessible from current data. The re- 
constructed estimations for this alternative second dark 
energy parameter are: 

{0.551 ± 2.079 for A/ = 0.4, 
-1.058 ± 0.835 for A/ = 0.6, (15) 
-0.733 ± 0.307 for A/ = 1.0 . 

This result shows more efficiently than the previous 
ones how difficult it is to extract information from small 



fitting windows. Moreover, the largest fitting window 
width we have chosen is again quite free of biasing ef- 
fects, and it gives values that are in agreement with those 
for the smallest window. The values corresponding to 
Ay = 0.6 seem to be consistent with a constant EoS, 
even though the errors on wq.s are really large. On the 
other hand, results from A^ = 1.0 are more compliant 
with a varying EoS. Both are also consistent with the re- 
sults in [23j, where the SNela Union2 sample was fitted 
with the {wo,wo,5) model. This is also a possible confir- 
mation that the reconstruction works fine, at least up to 
z w 0.4-0.5. 
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Figure 3. Second order polynomial reconstruction with SNela 
(Union2): reconstruction of the deceleration parameter q{z). 
DifTerent colors are used for different fitting window Af. from 
lightest to darkest for A/ = 0.4, A/ = 0.6 and A/ = 1.0 
respectively. The solid and the dashed black lines are respec- 
tively the WMAP7 ACDM (fi^ = 0.266) and the WMAP7 
(4- WMAP7.2) quiessence model. 
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^ Given the noisy behavior it would be more correct to speak about 
a first-time transition redshift, as it is possible that such line is 
crossed many times more because of the large uncertainties than 
a real physical source 

^ As it is well known, the CPL parametrization corresponds to 
w(z) =100-1- JiiaCl — a)- 



Figure 4. Second order polynomial reconstruction with SNela 
(Union2); reconstruction of the dark energy EoS w{z). Dif- 
ferent colors are used for different fitting window Af. from 
lightest to darkest for A/ = 0.4, A/ = 0.6 and A/ = 1.0 re- 
spectively. Blue contours correspond to the best fits for CPL 
from real data in 1281. 
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C. GRBs 

As discussed in the previous related section, the GRBs 
sample to be used is similar to the one studied in [6 , the 
difference being it is calibrated so that the GRBs Hub- 
ble diagram coincides with the SNela one in the redshift 
range where they overlap. 

We have performed a cut similar to the one we did 
for the SNela sample: the original data set was approx- 
imately homogeneous until a z sa 5 and shows only 2 
GRBs for higher redshift up to z sa 6.6. Trying to fit 
these data is pointless, and all one gets are large noisy 
trends at their related redshift; for this reason we perform 
fits without considering these last two points. 

The minimum length that can be used with GRBs for 
which the number of datapoints is higher than the num- 
ber of parameters is Aj = 1.3. Nevertheless, this is a 
border line case, because the number of datapoints is far 
too low and in many cases exactly equal to the number 
of fitting parameters, so that eventual results are statisti- 
cally invalid. This also gets clearly reflected in the error 
bars, which are rather big and give no information at all 
about the trends of the quantities. The minimal num- 
ber requirement we have also mentioned before (having 
> 10 points in each fitting window) is first fulfilled by 
Ay = 3.0. For further insight, we also present data from 
a wider fitting window with Aj = 4.0, while one must 
always keep in mind all the problems that enlarging A y 
way too much brings about. 

In Fig. [5] we can see the reconstructed variables y{z), 
y'{z), y"{z) for the GRBs sample. The reconstructed 
y{z) seems to be conditioned by the low number of data 
points in the sample, which translates into oscillations 
and growing error bars at high redshifts, larger than those 
for the SNela reconstruction. The first order derivative 
y' {z) presents two main differences with respect to the 
SNela case: 1. it decreases faster, and most importantly, 
2. its present value is sensibly lower. This will have a 
strong impact in the reconstruction of H{z), which has 
a strong dependence on y'{z) and becomes manifest in 
Fig. [sjd). This quantity displays a large deviation from 
current data, and it results in a very large Hubble con- 
stant value (with a large uncertainty as well): 



Hn = 



115 ±63 
126 ±51 



for 
for 



Af - 3.0, 
A/ = 4.0. 



(16) 



These values clearly differ from current accurate esti- 
mates from other astronomical probes '16], and no good 
agreement is found with ^18] as the global trend is con- 
cerned. 

When considering the second order derivative y"(z), 
represented in Fig. [5je) , one can notice the abrupt change 
at z « 0.8, which is also reflected in the related param- 
eters, i.e. q{z) and w{z), as shown in Figs. |5jf) and 
[5f^g). The deceleration and the EoS parameter trends 
are conditioned by very large error bars; the w profile, in 
particular, seems to be compatible with a constant EoS, 



w 1 for z < 0.8 and w - 0.1 for 0.8 < z < 2.; but 

also becomes divergent at z 3.5 due to the denominator 
entering its expression. 

The observed large discrepancy between the GRBs and 
the SNela reconstruction is probably both due to the 
low number of GRBs data in this regime and to intrinsic 
shortcoming in the calibration laws required for the treat- 
ment of GRBs datapoints. These disadvantage make us 
feel the GRBs are not reliable to be used in this model- 
independent framework. However, we expect that an im- 
provement on the knowledge of GRBs can probably lead 
to tighter measures and then to a better implementation 
in this method. 



D. Mock 

The main question arising from the previous section is: 
what is the degree of confidence on the reconstruction of 
cosmological quantities with this method? Put in other 
words, is there any way to discriminate and disentan- 
gle the "real" reconstruction from noisy features caused 
by shortcomings in the algorithm? And considering the 
particular kind of data we have used, this can also be 
translated into another question: until what redshift is 
the reconstruction well based? 

To try to throw some light into these matters, we have 
used a mock SNela sample: this is built on a fiducial cos- 
mological model, so that, in theory, if the method works 
well, one should expect to obtain reconstructed quantities 
compatible with this model. To create SNela mock sam- 
ples we have taken into account the specifications of the 
Wide-Field Infrared Survey Telescope (WFIRST) ^24^, a 
space mission planned by NASA, which has among its 
primary objectives to explore the nature of dark energy 
employing three distinct techniques: measurements of 
baryon acoustic oscillations, SNela distances, and weak 
gravitational lensing. The redshift distribution of the ex- 
pected SNela to be observed within various redshift bins 
in a survey of that sort is reported in |25j . The main red- 
shift range extends up to z ~ 1.6, with 4 expected SNela 
at z > 1.6. The very low redshift subsample, i.e. the 500 
SNela at z < 0.1, is supposed to come from the Nearby 
Supernova Factory (SNfactory) [2 6) with measurements 
in the redshift range (0.03,0.08), so we are going to as- 
sume the former as the lowest redshift value in our mock 
sample. 

Our formulae for errors on SNela magnitudes stems 
from a recipe used in the binning approach 127] , which we 
have however adapted for its application to non-binned 
data. According to this prescription errors are calculated 
as follows : 



syst 



(17) 



where: aint — 0.15 is the intrinsic dispersion in magni- 
tude per SNela, assumed to be constant and indepen- 
dent of redshift for all well-designed surveys; Cpec = 



10 




(b) 



(d) 



m. f. 



(e) 



(f) 



(g) 

Figure 5. Second order polynomial reconstruction with GRBs (Schaefer-Cardone). Coordinate distance y{z) and its derivatives 
{y'{z) and y"{z)), and the derived cosmological quantities: distance modulus i-i{z), Hubble function H{z), deceleration param- 
eter q(z) and the dark energy EoS w{z). Different colors are used for different fitting window A/: from lightest to darkest for 
Af — 3, Af = 4. respectively. 
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5(Ti,/ln(10)cz is the error due to the uncertainty in the 
SNela pecuUar velocity, with = 500 km/s, c the 
hght velocity and z the redshift for any SNela; asyst — 
0.02(1 + z)/(l + Zmax) is the floor uncertainty related to 
all the irreducible systematic errors with cannot be re- 
duced statistically by increasing the number of observa- 
tions. The value 0.02 comes from the fact that WFIRST 
is a space-based mission (it would be 0.05 for terrestrial- 
based observations); z^ax is the maximum observable 
redshift in the considered mission and this linear term 
in redshift is used to account for the dependence with 
redshift of many of the possible systematic error sources 
(for example the Malmquist bias or gravitational lens- 
ing effects). The fiducial model we have used is the 
wcdm-|-sz-|-lens case from WMAP7-year P], which has 
h = 0.75to ll, = 0.259l^;|5;|^, and is phantom-like 
with w = -l.Utoil- 

We will use a mock sample generated with that pre- 
scription to study the relation of the proposed recon- 
struction method with three different features: 1. data 
dispersion; 2. changes in the fitting window length and 
3. the range of validity of the method. 

The first point is addressed starting from the algorithm 
used to generate the mock sample. We take a set of val- 
ues of the cosmological parameters in the fiducial model 
which comprises these cases: the best fit, the best fit plus 
the 1(7 error, and the best fit minus the la error. We then 
construct all possible combinations of them, and then for 
these sets, and for any redshift in the sample, we find the 
largest possible variation in the distance modulus. We 
store this set of values as they are going to serve as our 
dispersion seeds. The next step is to build four mock 
samples as follows. We take the distance modulus corre- 
sponding to each redshift according to the fiducial model, 
and then we add to it a random value extracted from a 
Gaussian centered at zero and with a dispersion which is 
a number of times (0,0.25,0.5 or 0.75) the correspond- 
ing value from our dispersion seeds vector. The first case 
corresponds to mock data points that perfectly resem- 
ble the exact fiducial model, while he last value gives a 
more realistic dispersion similar to the from the present 
Union2 sample. 

For the sake of clarity, we will show results for a win- 
dow length Af = 0.6, which is the best case for the 
Union2 sample. In fact, it should also fare well with 
WFIRST (recall that with a sample with more points and 
more dense, the fitting window could be diminished with 
equally good final results). Moreover, we will restrict the 
discussion to the deceleration and EoS parameter, which 
are the two most sensitive parameters to eventual fea- 
tures. 

From Figs.[6]it is clear that the reconstruction method 
depends on the dispersion of the observed data. In fact, 
an effect of this sort should not come as a surprise, as 
there is always an intrinsic floor for this dispersion (de- 
pending on the physics of the observed phenomena and 
on the reduction procedure of observational data) . Con- 
sider, for instance the ideal case of null dispersion, as de- 



picted in Figs.|6ja) and|6jb), which means that the only 
element that makes the data differentiating from the fidu- 
cial model is the error on the distance modulus. For this 
case the reconstruction fails at low and high redshift. As 
discussed in previous section, at low redshift we have bi- 
ased estimations of the parameters due to the asymmetry 
of the fitting window, while at high redshift we have few 
points in each fitting window. As long as the dispersion 
grows, the deviation from the fiducial model grows too 
along with the noisy features, finally reaching the more 
realistic case (0.75). This one is shown in Figs.|6](g) and 
[6](h), where results are quite similar to those derived from 
Union2. It follows that an improvement on the observa- 
tional measures will be translated into a better recon- 
struction of the parameters under study. However, it has 
to be taken into account the limitation at high redshift 
where the density of points is lower. In Fig. [T] we can 
also verify how the fitting window infiuences results: in 
the ideal case of no dispersion in the data (Fig. [7|[a)), 
effectively the smaller the fitting window, the better the 
reconstruction, with only small deviation at the high red- 
shift tail, where the density of points drops. In the more 
realistic case in Fig.[7j^b), where a larger dispersion in the 
data is present, we can see how the choice of the fitting 
window length is crucial. A large window gives a small 
bias in the estimation of the parameters; on the other 
side, smaller windows give a better result, but the noise 
makes it difficult to understand what is the best choice 
and up to what redshift we have the best reconstruction. 



V. BINNED APPROACH 

For completeness, we have also carried out an inde- 
pendent analysis using a binned approach so as to get 
independent values of the magnitudes of interest at cer- 
tain values of the redshifts. This is, of course, opposed 
to the work spirit of the previous sections, which were 
devoted to the study of general trends of the H{z) or the 
dark energy properties q{z) or w{z). The procedure we 
present and discuss in this section can be used as an in- 
dependent check for theoretical models. We must stress, 
however, that for this part of our analysis we have con- 
sidered SNela only, given that according to our findings 
on GRBs they do not stand at the moment on the same 
quality grounds, at least in what concerns the reconstruc- 
tion methods considered in this paper. 

Our starting point has been to divide the SNela data 
sample into different numbers of redshift bins. Then, us- 
ing the data points within each of them, we have obtained 
y{z), y'{z) and y"{z). For this task we have followed ba- 
sically the steps described in the previous section, but 
with the slight modification of giving a low probability 
to the points near the boundaries of the bin. 

Results obtained with this approach are presented in 
Figs. |8] - [To] and in Table |v| For H{z) the binned ap- 
proach is well behaved, and the results agree with those 
presented in ^6 . However, for the dark energy EoS pa- 
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Figure 6. Deceleration q{z) and dark energy EoS w{z) - second order polynomial with WFIRST mock data with A/ = 0.6 and 
0.02. Black lines are the considered quantities in the fiducial cosmological model used to create the mock data set, the 



wcdm+sz+lens case from WMAP7-year [9], with h = 0.75to }4, ^2™ = 0.259lo:o9l and is phantom-hke with w = -1.12j;o;43. 
First panel: Null dispersion mock data. Second panel: Mock data dispersion factor 0.25. Third panel: Mock data dispersion 
factor 0.5. Fourth panel: Mock data dispersion factor 0.75. Dispersion values are as described in the related section. 
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N 


-2'med 


H{Zmed) 




qiZmed) 


7 bins 


77 
76 
77 
77 
76 
77 
76 


0.020 
0.037 
0.120 
0.262 
0.383 
0.554 
0.872 


72.376 ± 4.717 
89.058 ± 6.928 
74.221 ± 2.238 
84.428 ± 6.905 
72.376 ±4.717 
89.058 ± 6.928 
74.221 ± 2.238 


-21.507 ±59.869 
-17.654 ± 13.774 

-0.528 ± 2.985 
-10.709 ± 13.375 

-21.507 ±59.869 
-17.654 ± 13.774 
-0.528 ± 2.985 


-23.477 ±66.523 
-21.253 ± 16.817 

-0.036 ± 3.027 
-9.818 ±12.307 
-23.477 ± 66.523 
-21.253 ± 16.817 

-0.036 ± 3.027 


6 bins 


90 
89 
88 
88 
89 
90 


0.020 
0.043 
0.180 
0.327 
0.493 
0.817 


71.263 ±3.810 
77.258 ±2.966 
74.085 ± 3.027 
77.263 ± 6.982 
82.143 ±8.881 
108.020 ± 11.380 


-19.860 ±36.642 
-4.103 ±3.021 
-0.502 ±4.471 
-3.547 ± 15.760 
-3.861 ± 18.912 
—0.372 ± 5.025 


-21.380 ±40.009 
-4.170 ±3.312 

0.033 ±4.160 
-2.177 ± 11.795 
-1.664 ±10.408 

0.306 ± 2.510 


5 bins 


108 
107 
107 
106 
107 


0.021 
0.054 
0.254 
0.438 
0.769 


70.236 ±2.907 
70.792 ± 1.973 
81.099 ±3.568 
93.850 ± 9.979 
117.410 ± 12.170 


-13.304 ±23.791 
-0.528 ± 1.393 
-0.449 ± 5.014 

—3.413 ± 10.414 
-1.407 ±3.107 


-13.990 ± 25.630 
-0.057 ± 1.467 
0.083 ± 4.657 
—2.427 ± 8.840 
-0.535 ± 2.266 


4 bins 


134 
132 
132 
132 


0.023 
0.098 
0.357 
0.703 


69.465 ± 2.183 
71.411 ± 1.535 
82.531 ± 4.939 
109.330 ± 8.647 


5.196 ± 8.835 
-0.549 ± 0.976 

2.630 ± 5.835 
-1.077 ±2.158 


6.103 ± 9.345 
-0.052 ± 0.979 

2.607 ±4.287 
-0.269 ± 1.541 


3 bins 


177 
176 
176 


0.024 
0.216 
0.598 


72.199 ± 1.588 
75.981 ± 1.886 
97.006 ± 5.010 


-1.432 ± 1.553 
-0.980 ± 1.407 
-0.785 ± 1.288 


-1.086 ± 1.706 
-0.389 ± 1.274 
-0.030 ± 0.849 


2 bins 


264 
264 


0.026 
0.438 


71.913 ± 1.105 
89.085 ± 2.712 


-1.074 ±0.395 
-0.895 ± 0.496 


-0.684 ±0.418 
-0.204 ± 0.389 



Table I. Binned values for H{z), w{z) and q{z), where N and Zmed are respectively the number of data in each bin and the 
median redshift in each bin. 
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Figure 7. Dark energy EoS function 'w{z) - second order poly- 
nomial with WFIRST mock data. Left: mock data with null 
dispersion and different windows functions (from lightest to 
darkest, A/ = 0.4, A/ = 0.6 and A/ = 1.0). Right: mock 
data with dispersion 0.75 and different windows functions 
(from lightest to darkest, A/ = 0.4, A/ = 0.6 and A/ = 1.0). 



rameter 'w(z) and acceleration parameter q{z), we obtain 
large uncertainties when the number of bins is large. Nev- 
ertheless, these uncertainties get smaller as we decrease 
the number of independent bins and thus increase the 
amount data within a redshift bin. 



VI. SUMMARY AND CONCLUSIONS 

The abundant literature sources on model independent 
reconstructions of dark energy makes one believe this has 
become a worthy line of investigation. The approach we 
have been concerned with in this paper has been amply 
covered in [IHB] . The main idea is to take the dimension- 
less coordinate distance and try to reconstruct it with 
minimal assumptions, and then push the reconstruction 
two levels deeper and consider as well the first and sec- 
ond derivatives with respect to redshift. The way to do 
this is to carry out polynomial fits in different redshift 
windows (with carefully chosen characteristics). 

We have described the method following the original 
prescriptions, including throughout additional remarks 
to facilitate comprehension of the various steps involved. 
Our main work tool has also been SNela luminosity mea- 



surements, but we have been benefited by the availability 
of a sample which offers more than twice the number of 
sources: the very recent Union2 sample. This clearly al- 
lows for some improvement, as increased density of the 
datasets induces in general tighter constraints. Never- 
theless, it must be stressed that the sample considered in 
our paper is rather heterogeneous. 

A fact that seems to have been overlooked in previous 
references is that the quite noticeable oscillations arising 
at some levels of the reconstruction are mainly a conse- 
quence of data with large error bars in the samples. The 
fitting algorithm does not seem to be the culprit for these 
effects. We have been able to reach these conclusions by 
considering conveniently cut versions of the Union2 sam- 
ple. The cut produces a spectacular reduction of the 
noisy features. This suggests it could be worth revisit- 
ing this technique when significant advances occur in the 
compiled datasets. For completeness and comparison we 
have also performed a similar analysis with GRBs, but 
the reconstruction turn out to be rather poor (as other- 
wise expected). 

On general grounds we can say the reconstruction from 
SNela is very satisfactory for the dimensionless coordi- 
nate, but it degrades considerably as subsequent deriva- 
tives of it are considered. Estimates of global features 
agree very well with state of the art guesses from other 
sources, and no strong dependence on the size of the win- 
dow is appreciated at lowest levels. We find a Hubble 
constant of a bit more than 70 (km/s)/Mpc, a decelera- 
tion factor in the vicinity of —0.5, and a current value of 
the dark energy parameter in the vicinity of —1. 

For further insight into the problem we have also per- 
formed reconstructions using a mock SNela sample de- 
picting a future optimistic observational situation. The 
results support broadly what real data tell us, and in 
particular we find that dispersion in the data influences 
the reconstruction process a great deal. 

And flnally, using the binning method, we have ob- 
tained independent values of several quantities of inter- 
est. These may be useful for a quick/preliminary esti- 
mate of cosmological constraints without having to resort 
to a full dataset. 
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Figure 8. Values of H{z) for the binned approach for different values of independent bins, red points. The central point 
indicates the median redshift of the bin with its corresponding error bars (vertical lines). The horizontal bars delimite the 
redshift bin. The H{z) from 18 is also plotted (black points) in order to compare the results 
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